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ABSTRACT 

The T e ff location of Pre-Main Sequence (PMS) evolutionary tracks depends on the 
treatment of over-adiabaticity. We present here the PMS evolutionary tracks computed 
by using the mixing length theory of convection (MLT) in which the cvmlt = I / H p 
parameter calibration is based on 2D-hydrodynamical models (Ludwig et al. 1999). 
These MLT-a 2D stellar models and tracks are very similar to those computed with 
non-grey ATLAS9 atmospheric boundary conditions and Full Spectrum of Turbulence 
(FST) convection model both in the atmosphere and in the interior. The comparison 
of the new tracks with the location on the HR diagram of pre-MS binaries is not 
completely satisfactory, as some binary components are located at too low T e g. Be- 
sides, the pre-MS lithium depletion in the MLT- a 2D tracks is still much larger than 
that expected from the observations of lithium in young open clusters. This result is 
similar to that of FST models. Thus, in spite of the fact that 2D RHD models should 
provide a better convection description than any local model, their introduction is not 
sufficient to reconcile theory and observations. Lithium depletion in young clusters 
points towards a convection efficiency which, in pre-MS, should be smaller than in 
the MS. The pre-MS lithium depletion decreases significantly in FST models if we 
reduce the solar metallicity down to the value suggested by Asplund et al. (2004), but 
the corresponding solar model does not reproduce the depth of the convective zone as 
determined by helioseismology. 
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1 INTRODUCTION 

The location in the Hertzsprung-Russell diagram (HRD) 
of Pre-Main Sequence (PMS) evolutionary tracks is very 
sensitive to modeling details such as opacity, atmospheric 
boundary condit io ns, rota tion an d convection treatment. 
iMontalban et alJ (l2004al) and ID Antona fe Montalbanl 
i2003t) pointed out that, at least for T e ff down to 4000 K, 
the treatment of convection transport plays a role on the 
PMS evolutionary tracks that is much more important than 
the effects from recent improvements of low temperature 
opacities. They showed that it is the "average efficiency" of 
convection in the envelope which determines both, the T e a 
of PMS tracks and the lithium depletion 1 during this phase 
of stellar evolution. 



* E-mail: j.montalban@ulg.ac.be (JM); dantona@oa-roma.inaf.it 
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1 This conclusion is valid only if no other transport process dif- 
ferent than convection occurs during the PMS. 



Unfortunately, the approaches traditionally used in 
stellar evolution to describe the heat transport by con- 
vection are quite rough. The mixing length theory 
(MLT, Bohm-Vitense 1958) as well as the Full Spec- 
trum of Turbulence (FST) model llCanuto fc Mazzitehl 
Il99li ICanuto Goldman fc Mazzitellil Il99rj) are local theo- 
ries, which ignore physical processes such as convective over- 
shoot and radiative transfer effects. 

The efficiency of convection in the superadiabatic re- 
gion at the top of convective stellar envelopes determines 
the asymptotic value of the entropy in the deep and adia- 
batically stratified layers and, therefore, the star radius and 
T e g. Any treatment of convection, adjusted to obtain the 
correct adiabat, will provide similar global structures. This 
is the base of the MLT solar calibration. Given the rela- 
tionship between convective flux and over-adiabaticity in 
the MLT (-Fionv oc (V — Vad) 3 ^ 2 OMLT)i increasing the value 
of the omlt parameter we decrease the value of the over- 
adiabaticity, and the radius decreases. 

The only parameter in the MLT, lvmlt, is tuned to re- 
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produce the solar radius at the solar age. This method has 
provided very good results, in spite of its simplicity. Nev- 
ertheless, the Sun yields only one calibration value, and no 
physical principle guarantees this value to be appropriate 
for the whole HR diagram, or the whole stellar structure. 

The mass of the stellar convective envelope increases 
with decreasing T e g, and its top gradually shifts towards 
larger depths, so that the structure of the surface layers is 
less and less affected by convection as T e a decreases. For 
temperatures lower than 4700 K, however, the convection 
zone rises again due to H2 dissociation. Furthermore, as 
gravity decreases, the convective flux decreases because of 
the lower density, and therefore, the over-adiabatic region is 
more and more extended in the stellar atmosphere. Specific 
atmosphere models with an adequate treatment of radiation 
and convection transports are hence mandatory. 

The available grids of model atmospheres for general use 
of stellar structure (ATLAS9 by Kurucz 1993 and NextGen 
by Hauschildt et al. 1999) were generally computed just for 
one specific convection model, e.g. a given omlt = «mlt 
value in the MLT. Hence, in the solar calibration, only 
the MLT parameter use d in the interior can be changed. 
IMontalban et afl ]2004al) showed that this kind of procedure 
can lead to an uncertainty of the order of 200 K in the 1 Mq 
PMS track (without changing the effective temperature of 
MS), since its T c ff depends on the optical depth chosen to 
match the interior and the atmosphere models. On the other 
hand, Heiter et al. (2002) computed FST ATLAS9 models 
with a value of the «fst parameter fixed by a grey FST so- 
lar calibration. A non-grey FST solar calibration requires a 
larger value of q:fst (0.18 in the atmosphere and in the inte- 
rior -Ventura, private communication). Nevertheless, given 
the different properties of cvfst and cvmlt parameters, a 
discrepancy between a F ST and apST have no very serious 
consequences on the stellar structure o r on th e evolutionary 
tracks (as shown in Montal ban et al.l i2004al) . this discrep- 
ancy can lead to an uncertainty in T e g-ai given mass and 
luminosity- of the order of 2%, and <1% for 1 Mq). In any 
case, non-grey FST stellar models (either with cipsx^apg^, 
or with a P 1 ST =2-Q F t ^ r ) predict a too large pre-MS lithium 
depletion compared to observational data in open clusters. 

A more realistic approach consists in solving the hy- 
drodynamic equations coupled to the equation of radiative 
transfer. Recently much progress has been made in perform- 
ing 2-3D radiative-hydrodynamical (RHD) simulations of 
stellar surface convection (e.g. Freytag et al. 1996; Stein & 
Nordlund 1998; Asplund et al. 2000; Ludwig et al. 2002). 
These results, however, cannot be directly used in stellar 
evolution computations. A way of overcoming the problem 
of computing convection would be to have grids of 3D non lo- 
cal models, and, at each (T fr-gravity) , calibrate the a value 
which provides the same specific entropy jump between the 
atmosphere and the adiabatic region. This procedure does 
not give information on the structure of the overadiabatic 
layers, but allows a proper computation of the interior 2 . 

2 The use of an "average efficiency" of envelope convection is then 
useful if we wish to obtain general structural information con- 
cerning the stellar interior properties, exactly such as the lithium 
depletion, and the evolution of T c g and gravity. Of course, it is 
less useful if we need to know the structure of the overadiabatic 
layers, e.g. for the computation of acoustic oscillation modes. 



Extensive 3D RHD simulations for general applications 
to the computation of stellar evolution are not yet avail- 
able, but 2D RHD simulations have been performed by 
iLudwig et alJ £l999) for a large grid of T e g and gravity val- 
ues, and the structural information from these models has 
been translated into an effective mixing-length parameter 
q mlT) suitable to construct standard stellar structure mod- 
els. 

As the results from RHD models are in principle more 
reliable than the local models, we computed pre-MS and 
MS evolution, based on ILudwig et al] dl999li . for masses 
from 0.8 to 1.5 Mq. We compare the resulting evolutionary 
tra cks and Li abundance wi th those from recent local mod- 
els (IMontalban et ai]l2004aT) . which assume FST convection 
both in the non-grey atmosphere and in the envelope struc- 
ture. 

The application of 3D RHD atmosphere computation 
to studies of spectral line formation, have also shown that 
in many cases standard ID analyses are very misleading in 
terms of derived element abundances (Asplund 2005). In 
particular, 3D computations together with a much more 
accurate analysis of solar spectrum and a non-LTE treat- 
ment, led to a drastic revision of the photospheric solar 
composition. Asplund et al. (2004 and 2005) analysis re- 
sults in a significant reduction of solar metallicity down to 
Z/X ~ 0.017 - 0.018. We expect that a smaller lithium de- 
pletion could be obtained if we revise downward the solar 
metallicity. Unfortunately, the 2D RHD models by Ludwig 
et al. (1999), and the corresponding a 2 ^, calibration, are 
only available for the higher "old" solar metallicity and he- 
lium mass fraction Y=0.28. As we show in section 3 that 
FST and MLT- a 2D results are very similar for the standard 
"old" solar metallicity, we compute FST pre-MS tracks with 
the scaled down metal abundances, to test if the new abun- 
dances could solve the problem of PMS-lithium depletion in 
young open clusters. 



2 "AVERAGE EFFICIENCY" OF 

CONVECTION FOR PRE-MS MODELS 

ILudwig et al] (I1999T) . by using their 2D RHD atmosphere 
models, provided a calibration of the omlt parameter as a 
function of T c s and log g in the domain T c s — 4300 — 7100 K, 
logg = 2.54 — 4.74. We would like to stress two inter- 
esting aspects of this calibration in the PMS evolution 
context: first, the 2D models indicate that convection is 
on average 'efficient' in the atmosphere and the envelope, 
corresponding to a large «mlt value; and, second, Qmlt 
increases as T c g decreases. The behavior has been 

confirmed by 3D atmosphere models. Ludwig, Allard and 
Hauschildt (2002) find, for an M dwarf at T cff = 2800 K and 
logg = 5, cvmlt— 2.1, a value n oticeably larger than that 
required for the Sun faMT,T~1.6l. iTrampedach et al] (ll99STl 
made a similar a® 1T calibration by using 3D model atmo- 
spheres computed by means of a different code (the time- 
depend ent, compressible, explici t, radiative- hydrodynamics 
code bv lStein fc Nordlundlll99cfl . For the range of main se- 
quence gravities and log 2^1=3.68-3.83, this cimlt calibra- 
tion shows the same behavior than the a™ LT one and a 
systematic offset, in the sense a^ LT > o^ LT . For the Sun, 
ILudwig et al] ]l999h found than their 2D-based calibration 
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Figure 1. Evolutionary tracks using FST in the logT e ff vs. log g 
plane (solid line; non— grey models with T p h = 10 by Montalban 
et al. (2004a) and 2D calibrated MLT (dashed line). 



underestimates the value of qmlt by a factor ~ 0.2, and 
they explained that as due to a combined effect of low- 
temperature op_a£ities_and t he 2D approximation. On the 
other hand, lAsolund et alJ ((2(3)00) have compared 2D and 
3D atmosphere models for Sun, and find that the 2D so- 
lar model has margi nally larger gradients than the 3D one. 
iLudwig et alJ (l999) propose the use of a constant scaling 
factor to compensate the systematic offset seen in the Sun, 
as the relative variations of qmlt with T e s and log g should 
be less affected by systematic shortcomings such as the 2D 
approximation. 

It has also been shown that there is no unique value 
of the qmlt parameter that can reproduc e the temperature 
gradi e nts in the ove r -adiab atic regions dSteffen fc Ludwiel 
Il999l : iTrampedachl 12004). Nev e rthele ss, the function 
Qmlt (T e g, log g) by ILudwig et all lll999l) . guarantees that 
the adiabat of a stellar model, computed with this a™ LT 
and grey boundary conditions, provides the same adiabat 
(and therefore the same radius and effective temperature) 
as one from a complete 2D RHD computation. 

An additional advantage of the procedure of comput- 
ing MLT-q 2D models is that the entropy in the adiabatic 
convection region is independently determined and is not af- 
fected by the ambiguity and the uncertainties deriving from 
the use of several different physics in the over-adiabatic re- 
gion of the star. So, standard non-grey MLT models are not 
necessarily equivalent to MLT-2D/3D ones even if the nu- 
merical value of lvmlt parameter used in the stellar structure 
computations is in both cases the same or of the same order. 
We must indeed be careful when looking at the exact mean- 
ing of published results. For instance, the MLT- a 2D models 
will not be similar to t he popular set of no n-grey models 
having "qmlt = 1.9" bv lBaraffe et alJ <ll998l) . In fact, these 
models have indeed qmlt = Qmlt = 1-9 in the computa- 
tion of the internal structure, but the atmosphere down to 



T p h = 100, i.e. to the match point with with the interior, 
is computed with cvmlt = 1. It has been ofte n pointed out 
that the PMS tracks bv lBaraffe et alJ il99ST) match better 
observations since they have lower T e ff than other compu- 
tations (e.g. Hillenbrand & White 2004) and, sometimes, 
that has been attributed to the improved opacity in the 
adopted NextGen atmospheres (Allard & Hauschidt, 1997, 
hereinafter AH97). On the contrary, for mass values close to 
the solar mass, a large part of this result is simply due to the 
very low efficiency of convection in the most superadiabatic 
part of the star, namely the atmosphere down to r p h = 100 
(Montalban et al. 2004a; D'Antona fc Mont alban 2003). 
The calibration of Qmlt in ILudwig et alJ lll999l) . whose zero 
point must be obtained by reproducing the solar radius, is 
valid for grey models, so the same Qmlt value is used in 
the whole structure, and the results are less ambiguous and 
easier to be interpreted. 



3 MLT— a MODELS 

The MLT-q 2D stellar models have been computed with 
the code ATON2.0 (Ventura et al. 1998a), with grey 
boundary condition s, and the MLT option for convection. 
ILudwig et alJ dl999l) suggest, in the context of stellar evolu- 
tionary models, to calibrate Qmlt with the present Sun, and 
to use its ratio to Qmlt(Q) as scaling factor for the func- 
tion Q M ^ T (T e ff,log g). So, we used the value Qmlt — 1-6 (ob- 
tained from the grey MLT Sun calibration wit h ATON2.0) to 
scale t he analytical fits of qmlt (T e g , log g) bv ILudwig et alJ 
(1999). These "solar-calibrated" analytical fits were intro- 
duced in ATON2.0 to allow for a continuous variation of the 
Qmlt parameter as the star evolves in the HR diagram. 

The adopted helium mass fraction is Y=0.28 and the 
metal mass fractio n is Z=0.02. The other physical inputs 
are the same as in lMontalban et alJ <|2004af) 3 . In Fig. Qwe 
superimpose the MLT-q 2 D tracks (dashed - lines) and the 
non-grey FST models bv iMontalban et al] J2004all in the 
pla ne log g vs. T e g. Th e latter tracks are computed by using 
the lHeiter et al] <l2002l) boundary conditions, and adopt FST 
convection, both in the atmosphere and in the interior. The 
two sets of results, in the range 0.8 - 1.5M® are very similar. 
Since stellar radii (and therefore the Tiff's) are determined 
by the efficiency of convection in the over-adiabatic region, 
we can conclude that the global efficiency of convection in 
MLT- a 2D and FST models is very similar. This allows the 
use of FST models as proxies of 2D RHD models, also for 
different chemistry. So, we have computed low metallicity 
FST models with chemical and convection parameters given 
by the new solar calibration, that is: Y=0.2495, Z—0. 01305 
and the FST qfst parameter equal to 0.117. 

3.1 Pre-MS Lithium depletion in MLT-q 2D 
models 

Lithium is a light element that is burned by the nu- 
clear reaction 7 Li(p, 4 He) 4 He at a relatively low tempera- 
ture (T Li ~ 2.5 10 6 K). The PMS lithium depletion rate 

3 The solar mixture used in the opacity and equation of state 
tables is that from Grevesse &: Noels 1993 
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Figure 2. Lithium evolution for models computed with 2Da cal- 
ibrated MLT (solid lines) and complete (non grey) FST models 
(dashed lines) 




log t 



Figure 3. Lithium evolution for the solar mass with different as- 
sumptions about convection and model atmospheres. The dotted 
line at bottom represents today's solar lithium abundance. MLT 
models with AH97 model atmospheres down to r p h = 10 and 100 
are shown dotted for oj; n = 1 and dash- dotted for a; n = 1.9. The 
Montalban et al. (2004a) MLT models with Heiter et al. (2002) at- 
mospheres down to r p h = 10 (lower) and 100 (upper) are dashed; 
The continuous lines show the non-grey FST models for r p h = 10 
and 100, and, in between, the long dashed line shows the MLT- 
o 2D model. 



strongly depends on temperature and density at the bot- 
tom of the convective region. Lithium abundances detected 
in young open clusters , such as the Pleiades and cvPer (e.g. 
ISoderblom et alll99&t) . indicate that solar mass stars do not 
deplete a significant fraction lithium during their PMS. That 
implies structures of these stars with a temperature at the 
base of their convective envelope smaller than TLi. 

It is generally accepted (Gough fc Weiss 1976; 
Christensen-Darlsgaard 1997) that the overall late-type 
stars structure does not depend on the details of the treat- 
ment of over-adiabatic regions. Hence, the radius is de- 
termined by the value of the entropy jump between the 
photosphere and the deep adiabatic regions. Christensen- 
Dalsgaard (1997) pointed out that, for the Sun, the depth 
of the convective zone is almost insensitive to changes in 
the adiabat of the convective zone, and that the changes in 
the surface radius are due only to the changes of the radius 
of the bottom of the convective zone. He also pointed out 
that the same should probably occur for other stars, if the 
derivatives of opacity with respect to temperature and den- 
sity are of the same order than in the Sun. The similarity 
between FST and MLT- a 2D -track locations suggests that 
the "average efficiency" of convection is quite close in both 
series of models, and hence we expect that also the lithium 
depletion in the MLT- cv 2D models should be similar to that 
of the FST ones. In Fig. [5] we plot the lithium depletion as 
a function of time for masses between 0.8 and 1.2 Mq. In 
these computations the initial lithium abundance was taken 
as log N(Li)=3.31, i.e . the s olar system abundance given 
bv lAnders fc Grevessei (Il989l) 4 . Both the MLT- a 2D mod- 
els and the non-grey FST models show a depletion of the 
order of ~ 1.7 dex for the 1 Mq evolutionary track. Thus 
the non-grey FST models provide a description of the stel- 
lar structure very similar to the RHD 2D models, also with 
respect to lithium depletion. 

Generally, local convection models are calibrated by 
requiring that the free parameter (s) reproduces the solar 
radius at the solar age. The MLT models satisfying this 
constraint predict, however, a pre-MS lithium depletion 
which is not compatible with that obse rved in young open 
clusters dP'Antona fc Montalblr] I2003D . In Fig. El we plot 
the lithium vs. t ime ev olution in some IMq models by 
Montalban et alJ J2004af) and in the MLT- a 2D models. As 
D'Antona fc MontalbanH2005t) remarked, the models adopt- 
ing AH97 atmospheres and Qmlt — 1 deplete lithium in pre- 
MS by only 0.15 dex. The same models with Q!ml T =1.9 de- 
plete 0.7 dex if the matching point is r p h=100, but more 
than 1 dex for r p h=10. We stress again that this result is 
only due to the use of a less efficient convection when the 
match with the atmosphere (^mlt = ^) is done at r p h=100, 
even if the Qmlt parameter chosen for the interior, in or- 
der to fit the solar radiu s, is much larger (fo r instance 1.9). 
As it was shown in lMontalban et al . (2004a) and recalled in 
Sect. 1, the choice of the optical depth at which we match 
the two parts of the stellar structure, is not without conse- 
quences. It implies an uncertainty of almost 200 K for the 



4 The curre nt adopted value for the meteoritic Li abundance is 
3.25 ± 0.06 iAsplund et alj|2 005b1. but the small difference with 
respect to Anders &: Grevesse (1990) does not significantly affect 
our results 
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T e ff of 1 Mq PMS. The optical depth at which the external 
layers change from adiabatic to over-adiabatic stratification 
depends on T e a, logg, and chemical composition as well. To 
match an qmlt = 1 structure with an ayiur=l-9 one at 
fixed T p h along the whole evolution is equivalent to decrease 
the average efficiency of convection in a complicate, hidden 
and unjustified way. 

Li abundance in young clusters suggests that, if no other 
physical process is affecting the location of the base of the 
convective zone, the PMS-tracks should be cooler than solar 
calibrated convection predicts. This evidence should be used 
to look for the shortcomings of our stellar models. 

3.2 Pre-MS binaries 

The location of pre-MS binaries in the HR diagram is a 
powerful way of constraining stellar models, but unfortu- 
nately its indications are, by now, still ambiguous, because 
of the dearth of systems with well known masses and and 
good atmospheric parameters. For some of these binaries, 
the location is consistent with models having high convec- 
tion efficiency, but for others this is not the case. In Figure^] 
we compare the MLT-« 2D tracks with the parameters de- 
rived from the observations of four binary systems: RXJ 
05 29.4+0041 (accordin g to the most recent determination 
bv ICovino et alJ J2004M . V1174 Ori (Stassun et al. 2004), 
NTT 045251+3016 (Steffen et al. 2001) and HD 98800 B 
feoden et aljl2005h . A detailed analysis is out of the scope 
of this paper, but it is evident that the agreement is reason- 
ably good for most of the primaries but much less satisfac- 
tory for the secondary, i.e. less massive, components. The 
VI 174 Ori primary mass (IMq) is consistent with the cor- 
responding MLT-« 2D track, as well as the primary masses 
of RXJ 0529.4+0041 (1.3M ) and HD 98800 B (0.699M Q ), 
but that of NTT 045251+3016 is rather different from the 
expected value (1.45+0. 19Mq ) is 5 . Besides, at least two 
secondaries seem to be too cool for the MLT-a 2D tracks, 
and one could perhaps guess a systematic behavior, in the 
sense that the low mass PMS stars seem to be cooler and to 
have a larger than expected radius. 

This recalls what happens for the few known late-type 
low-mass MS binary components. High quality observations 
(eclipsing binaries and/or interferometric measurements) re- 
vealed a meaningful discrepancy between stellar parameters 
derived from observations and those based on theoretical 
models: the theoretical stellar radii seem to be be under- 
estimated by 10-20%, and the effective temperatures over- 
estimated by 5% (200 K) (e.g. Ribas 2006; Torres et al. 
2005; and references therein). An accepted, but not defi- 
nitely proved, explanation is that these larger than expected 
stellar radii are related to stellar activity. In fact, the mag- 
netic fields associated with stellar activity can decrease the 
efficiency of convection (Gough & Tayler 1966; Stein et al. 
1992 and references therein) and as a consequence the stellar 
radius must grow to transport the same quantity of energy. 
Besides, the decrease in effective temperature could be due 

5 the same Stcf fen et alJ i200ll) remark that the location of the 
prima ry in this system is 3cr away from the lD'Antona &: Mazzitellil 
jl994l) grey-FST tracks, which are hotter than the MLT-a 2D 
tracks by ~270 K. 
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Figure 4. The location in the HR diagram of the MLT-a 
tracks is shown together with four binaries with well determined 
masses, labeled in th e figure. Full squares (red): RX J 0529.4+0041 
iCovino et al.l2004l) : full circles (blue): V1174 Ori JStassun et alJ 
1 2004) : open squares (black): NTT 045251+ 3016 iSteffen et all 
boo II) : full triangles (magenta): HD 98800 B feoden et al.ll2005l) 

to spots on the stellar surface. In a very recent paper, Tor- 
res et al. (2005) have shown that the discrepancy between 
theory and observation can also appear for stars quite simi- 
lar to the Sun. In their study of the eclipsing binary V1061 
Cygni with M Aa = 1.282M and M Ab = 0.9315M©, they 
conclude: "Current stellar evolution models that use a mix- 
ing length parameter olml appropriate for the Sun agree well 
with the properties of the primary, but show a very large di- 
vergence in the radius of the secondary, in the sense that 
the predicted values are ~ 10% smaller than observed. In 
addition, the temperature is cooler than predicted by some 
200 K ... for a star only 7% less massive than the Sun". 
They suggest that there must be a relation among activity 
level, decrease of T c £f and underestimation of the stellar radii 
by the standard stellar models (calibrated by matching the 
Sun). 



4 FST MODELS WITH REVISED SOLAR 
METALLICITY 

We have considered, until now, MLT-« 2D models computed 
for the "standard" solar metallicity Z=0.02, but what will 
happen if we decrease Z to ^0.0 13, as suggested by t he re- 
analysis of the solar spectrum bv lAsplund et al. (2004)? The 
lithium depletion would probably be reduced, without af- 
fecting too much the pre-MS tracks location, as we still 
must calibrate convection in order to reproduce the solar 
location. RHD 2D models are not available for this revised 
metallicity, nevertheless, we can compute models with FST 
convection, and make the hypothesis that, for the reduced 
solar metallicity, they provide similar results to MLT-a 2D 
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models, as they do for the standard solar metallicity. In fact, 
the effective temperature in the Hayashi track corresponding 
to the solar model, calibrated with the new solar metallicity, 
is only 0.3% lower than that for the standard solar metal- 
licity. However, the temperature and density at the bottom 
of the convective zone of these PMS models, have changed 
by 6 and 12 % respectively. As a consequence, lithium has 
been depleted only by 0.5 dex. We recall, however, that the 
solar model with new solar metallicity is not able to fit ei- 
ther the depth of the convective zone, or the sound speed 
profile inside the Sun. 

Since the new solar abundances were published, several 
teams tried to recover the good match between the helioseis- 
mic Sun and the standard solar model. They proposed to 
increase the opacity, the microscopic diffusion coefficients, 
or a combination of both (e.g. Basu & Antia 2004, Mon- 
talban et al. 2004b, Guzik et al. 2005). We recall, that none 
of the proposed solutions is fully satisfactory. Nevertheless, 
for what concerns the lithium problem, whatever solution be 
adopted, it should be able to change the bottom of the con- 
vective zone during MS without changing the PMS models. 
If the discrepancy between the standard solar model and 
helioseismology could be solved only changing the micro- 
scopic diffusion (not reliable at the moment, see for instance 
Guzik et al. 2005), since this process is very slow, perhaps 
we could keep a low lithium depletion. For other processes, 
we should justify why they work during MS and not during 
solar PMS. Basu & Antia (2005) and Bahcall et al. (2005) 
proposed the solution of increasing the Ne abundance (not 
directly observed in the solar spectrum) by a factor ~ 3.5. 
This suggestion has been supported by Drake & Testa (2005) 
which, on the basis of Ne abundance determination in ac- 
tive stars, argue that the solar Ne could be 2.5 times larger 
than the value adopted for the solar mixtures. Without en- 
tering in the controvers y about these new Ne abundance 
determination (see e .g. ISchmelz et al. 1 120051 lYounel l2005t 
lAsplund et"al]|2005d) . we would like to recall that Ne is, 
with O, one of the main contributors to opacity at the bot- 
tom of the solar convective zone (close to TLi). The solution 
proposed by Basu & Antia (2005) implies to increase the 
opacity by replacing the decrease of O abundance in the 
new solar mixture (Asplund et al. 2004) by the increase of 
Ne one. A direct comparison between OPAL opacity tables, 
with GN93 mixture, and a solar mixture with Ne increased 
by a factor 3.5 yields that the opacity from this latter ta- 
ble, for the temperature and density typical of PMS stars, 
is even larger than for GN93 mixture. As a consequence, 
one can predict a deeper convective region in the PMS and, 
therefore, a lithium depletion even larger than in the "old" 
standard solar model. 



5 CONCLUSIONS 

The computation of MLT pre-MS tracks in which the omlt 
parameter is calibrated on 2D RHD models (MLT-a 2D 
tracks) shows that these tracks deplete too much lithium, 
at variance with the observations of stars in young open 
clusters. The MLT-q 2D tracks are similar, in HR diagram 
loc ation and lithium depleti on, to the FST non-grey tracks 
bv iMontalban et ail (|2004al ). Thus we compute FST non- 
grey models, correponding to the new solar metallicity de- 



termined bv lAsplund et ail J2004), as proxies of possible 2D 
RHD models for the reduced solar metallicity. The lithium 
depletion in 1 Mq model has been reduced from 1.7 dex 
(with the "old" solar metallicity) to ~ 0.5 dex. We know, 
however, that such a low metalllicity is not able to reproduce 
other, well established, properties of the Sun. Hence, much 
more work shall be done, analyzing all the consequences of 
the new solar metallicity, before adopting the new abun- 
dances as solution of the Li problem. 

The problem of lithium is therefore not solved in the 
framework of the numerical simulations of convection. On 
the contrary, 2D and 3D numerical simulations imply even 
more efficient convection in the PMS, while lithium deple- 
tion in young clusters and PMS binaries studied here suggest 
a lower efficiency. We conclude that convection in pre-MS 
must be less efficient than what is suggested by 2D RHD 
models. It is still possible that the description of pre-MS 
convection requires the introduction of a second parame- 
ter -linked to the stellar rotation and magnetic field, as we 
have suggested in the past (Ventura et al. 1998b; D'Antona 
et al. 2000). 

The data on masses and radii coming from high quality 
observations of eclipsing binaries and/or interferometry of 
late type stars agree only in part with the MLT-a 2D models. 
So, while some observations (stellar parameters of binaries 
and lithium abundances in young clusters) require a low 
efficiency of convection in the late-type domain of the HR 
diagram, numerical simulations of convection (Ludwig et al. 
1999; Trampedach et al. 1999; Ludwig et al. 2000) show the 
opposite: an equivalent value of the qmlt parameter that 
increases as effective temperature decreases. 

The abovementioned observational results, together 
with the high activity level that PMS-stars can show (see 
e.g. Tayler 1987) , seems to support the suggestion that that 
other physical processes (rotation, magnetic fields) affect the 
efficiency of convection in late type stars. 
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